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Abstract — The aim of this paper is twofold. On one hand, it 
presents a multi-dimensional Markovian state transition model 
characterizing the behavior at the Medium Access Control 
(MAC) layer by including transmission states that account for 
packet transmission failures due to errors caused by propagation 
through the channel, along with a state characterizing the system 
when there are no packets to be transmitted in the queue of 
a station (to model non-saturated traffic conditions). On the 
other hand, it provides a throughput analysis of the IEEE 
802.11 protocol at the data link layer in both saturated and 
non-saturated traffic conditions taking into account the impact 
of both transmission channel and multirate transmission in 
Rayleigh fading environment. Simulation results closely match 
the theoretical derivations confirming the effectiveness of the 
proposed model. 



I. Introduction 

The IEEE802.il MAC [1] presents a mandatory op- 
tion, namely the Distributed Coordination Function (DCF), a 
medium access mechanism based on the CSMA/CA access 
method, that has received considerably attention in the past 
years [2]-[14]. 

A number of papers [3]-[5], after the seminal work by 
Bianchi, have addressed the problem of modelling the DCF in 
a variety of traffic load and channel transmission conditions. 
Most of them focuses on a scenario presenting N saturated 
stations that transmit toward a common access point (AP) 
under the hypotheses that the packet rates along with the 
probability of transmission in a randomly chosen slot time 
is common to all the involved stations, while error events on 
the transmitted packets are mainly due to collisions between 
packets belonging to different stations. 

Modeling of the DCF of IEEE 802.11 WLANs in unsat- 
urated traffic conditions has been analyzed in a number of 
papers [6]-[10]. In [6] the authors extended the underlying 
model in order to consider non-saturated traffic conditions by 
introducing a new state, not present in the original Bianchi's 
model, accounting for the case in which the station queue is 
empty after successful completion of a packet transmission. 
Paper [7] proposes an extension of the Bianchi's model 
considering a new state for each backoff stage accounting for 
the absence of new packets to be transmitted, i.e., in unloaded 
traffic conditions. 
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In [11], the authors look at the impact of channel induced 
errors and the received SNR on the achievable throughput 
in a system with rate adaptation whereby the transmission 
rate of the terminal is adapted based on either direct or 
indirect measurements of the link quality. In [12], authors 
observed that in multirate networks the aggregate throughput 
is strongly influenced by the bit rate of the slowest contending 
station: such a phenomenon is termed performance anomaly 
of the DCF of the IEEE 802.11 protocol. In [14], authors 
provide DCF models for finite load sources with multirate 
capabilities, while in [13] authors propose a DCF model for 
multirate networks and derive the saturation throughput. In 
both previous works, packet errors are only due to collisions 
between different contending stations. 

In this paper, we substantially extend a previous work 
proposed in the companion papers [8] -[9] considering real 
channel conditions, both saturated and non-saturated traffic, 
and multirate capabilities. As a reference standard, we use 
network parameters belonging to the IEEE802.11b protocol, 
even though the proposed mathematical model holds for any 
flavor of the IEEE802.il family or other wireless protocols 
with similar MAC layer functionality. 

This paper is organized as follows. After a brief review 
of the functionalities of the contention window procedure 
at MAC layer, section [TT] substantially extends the Markov 
model initially proposed by Bianchi, presenting modifications 
that account for transmission errors. Section [ill] provides an 
expression for the aggregate throughput of the link, while 
Section[IV]derives the time slot duration needed for throughput 
evaluation. The adopted traffic model is discussed in Sec- 
tion [V] Section [VI] briefly addresses the modelling of the 
physical layer of IEEE 802.11b in a variety of channel fading 
conditions. In section \VU\ we present simulation results where 
typical MAC layer parameters for IEEE802.11b are used to 
obtain the throughput as a function of various system level 
parameters, and the SNR under typical traffic conditions. 

II. Markovian Model 

In a previous paper [8], we proposed a bi-dimensional 
Markov model for characterizing the behavior of the DCF 
under a variety of real traffic conditions, both non-saturated 
and saturated traffic load, with packet queues of small sizes, 
and considered the IEEE 802.11b protocol with the basic 2- 
way handshaking mechanism. Many of the basic hypotheses 
are the same as the ones adopted by Bianchi in the seminal 
paper [2]. 



As a starting point for the derivations which follow, we 
adopt the bi-dimensional model proposed in the companion 
paper [8], appropriately modified in order to account for a 
scenario of N contending stations each one employing a 
specific bit rate and a different transmission packet rate. For 
conciseness, we invite the interested reader to refer to [8]- 
[9] for many details on the considered bi-dimensional Markov 
model. 

Consider the following scenario: N stations transmit toward 
a common AP, whereby each station, characterized by an own 
traffic load, can access the channel using a data rate in the set 
{1,2, 5.5, 11} Mbps depending on channel conditions. Any bit 
rate is associated with a different modulation format, whereas 
the basic rate is 1 Mbps with DBPSK modulation (2Mbps with 
DQPSK if short preamble is used) [1]. We identify a generic 
station with the index s G S = {1, 2, • • • , N}, where N is the 
number of stations in the network, and S is the set of station 
indexes. As far as the transmission data rate is concerned, we 
define four rate-classes identified by a rate-class identifier r 
taking values in the set R = {1, 2, 3, 4} ordered by increasing 
data rates Rn = {1,2,5.5,11} Mbps (as an example, rate- 
class r = 3 is related to the bit rate 5.5 Mbps). Concerning 
control packets and PLCP header transmissions, the basic rate 
is identified by Rc. 

The traffic load of the s-th station is identified by a packet 
arrival rate (PAR) evaluated in packets per second. Upon 
defining both rate-classes and traffic, we can associate a 
generic station s with a rate-class r G R and a proper traffic 
load \^ s \ Therefore, we need to specify, with respect to 
the model proposed in [8], specific probabilities along with 
different Markov chains for each contending station in the 
network. 

The two sources of errors on the transmitted packets are 
both collisions between packets and channel induced errors. In 
relation to the s-th station in the network, collisions can occur 

(s) 

with probability P co /, while transmission errors due to imper- 

(s) (s) 

feet channel can occur with probability P e v . Notice that P e 
depends upon the station rate-class r, which is in turn related to 
the received Signal-To-Noise (SNR) (appropriate expressions 
will be provided in Section|Vl]for each rate-class). We assume 
that collisions and transmission error events are statistically 
independent. In this scenario, a packet from the s-th station 
is successfully transmitted if there is no collision (this event 
has probability 1 — P^fi ) and the packet encounters no channel 

(s) 

errors during transmission (this event has probability 1 — P e ). 
The probability of successful transmission is therefore equal 
to (1 — Pe )(1 — P^j), while the equivalent probability of 
failed transmission is defined as 

p(s) _ p(<0 . p(s) _ p(s) . pO) m 

r e-q ~ r col ^ r e r e e col w 

To simplify the analysis, we make the assumption that the 
impact of channel induced errors on the packet headers are 
negligible because of their short length with respect to the 
data pay load size [8]. 

The modified Markov model related to the s-th contending 
station is depicted in Fig. Q] We consider (m + 1) different 
backoff stages including the zero-th stage. The maximum 
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Fig. 1. Markov chain for the contention model of a generic station s 
in general traffic conditions, based on the 2-way handshaking technique, 
considering the effects of channel induced errors. 

contention windows (CW) size is W max — 2 m Wo, and the 
notation Wi — 2 1 Wq is used to define the i th contention 
window size (Wo is the minimum contention window size). 
A packet transmission is attempted only in the (t, 0) states, 
Vi = 0, . . . , to. In case of collision, or due to the fact that 
transmission is unsuccessful because of channel errors, the 
backoff stage is incremented, so that the new state can be 
k) with uniform probability P e g /Wi+i. The contention 
window is supposed to be common to all the stations in the 
network; for this reason the station index s is dropped from 
the contention model depicted in Fig. [T] 

In order to account for non-saturated traffic conditions, 
we introduced a new state labelled /, for the following two 
situations: 

> Immediately after a successful transmission, the queue 
of the transmitting station is empty. This event occurs 
with probability (1 — <z^)(l — Pe<i), whereby is the 
probability that there is at least one packet in the queue 
after a successful transmission, 
• The station is in the idle state with an empty queue 
until a new packet arrival in the queue. Probability pj S Q 
represents the probability that while the station resides in 
the idle state / there is at least one packet arrival, and a 
new backoff procedure is scheduled. 

We notice that the probability Pj o of residing in the idle state 

is strictly related to the adopted traffic model. 

III. Markovian Process Analysis 

This section focuses on the evaluation of the stationary state 
distribution of the Markov model proposed in the previous 
section. The objective is to find the probability that a station 
occupies a given state at any discrete time slot along with 
the stationary probability bj of being in the idle state. This 
mathematical derivation is at the basis for the derivation of 
the probability t^* 1 that a station will attempt transmission 
in a randomly chosen slot time. For the sake of simplifying 



the notation, in what follows we will omit the apex (s) since 
the mathematical derivations are valid for any contending 
station s = 1, . . . , N. For future developments, from the model 
depicted in Fig. [T] we note the following relations: 
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whereby bij is the stationary probability to be in the state 
labelled i,j of the Markov chain in Fig.Q] 

Let us focus on the meaning of the idle state / noted in 
Fig. [T| to which the stationary probability bj is associated. 
It considers both the situation in which after a successful 
transmission there are no packets to be transmitted in the 
station queue, and the situation in which the packet queue is 
empty and the station is waiting for new packet arrivals. The 
stationary probability of being in state bj can be evaluated as 
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Upon employing the probabilities 6,.o noted in (0, it is 
straightforward to obtain: 
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By using the previous result, (0 simplifies to 

, 1-9 



^0,0 



(4) 



(5) 



1,0 



The other stationary probabilities for any k £ [1, Wj— 1] follow 
by resorting to the state transition diagram shown in Fig. [T] 

q(l-P eq ) -J2Za ho+ 
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Employing the normalization condition, after some mathe- 
matical manipulations, and remembering @, it is possible to 
obtain: 
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From (0, the following equation for computation of 60,0 easily 
follows: 

bo.o = l -± (9) 



Equ. (0 is used to compute t^ s \ the probability that the s-th 
station starts a transmission in a randomly chosen time slot. 
In fact, taking into account that a packet transmission occurs 
when the backoff counter reaches zero, we have: 
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whereby we re-introduced the apex (s) since this expression 
will be used in the following. 

The collision probability needed to compute 

can be found considering that using a 2-way hand-shaking 
mechanism, a packet from a transmitting station encounters a 
collision if in a given time slot, at least one of the remaining 
(N — 1) stations transmits simultaneously one packet. Since 
each station has its own t^ s \ the collision probability for the s- 
th contending station depends on the transmission probabilities 
of the remaining stations as follows: 
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Given the set of TV equations ([TJ and ( TTOt . a non-linear system 
of 2N equations can be solved in order to determine the 
values of and P^J for any s — l,...,N: this is the 
operating point corresponding to the N stations in the network, 
needed in order to determine the aggregate throughput of the 
network, defined as the fraction of time the channel is used to 
successfully transmit payload bits: 



N 1 
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(12) 



whereby the summation is over the throughput related to the 
N contending stations, PL is the average payload size, and 
T av is the expected time per slot defined in the following. 

Probabilities involved in (fT2l i are as follows: P^ is the 
PER (or FER) of the s-th station due to imperfect channel 

(si 

transmissions, and P s v is the probability that a packet trans- 
mission from the s-th station is successful. In the next section, 
we derive the mathematical relations defining both T av and the 
probabilities involved in ( fT2l . 

IV. Estimating the Average Time Slot Duration 

In order to proceed further, we need to evaluate the average 
time T av spent by a station in any possible state. The average 
duration T av of a time slot can be evaluated by weighting 
the times spent by a station in a particular state with the 
probability of being in that state. It is possible to note four 
kind of time slots. The average idle slot duration, identified 
by Tj, in which no station is transmitting over the channel. 
The average collision slot duration, identified by Tc, in which 
more than one station is attempting to access the channel. The 
average duration of the slot due to erroneous transmissions 
because of imperfect channel conditions, identified by Tg. 
The average slot duration of a successful transmission, 
identified by T$. 

The average idle slot duration. The average idle slot 
duration can be evaluated as the probability (1 — P t ) that no 
station is attempting to gain the access to the channel times 
the duration a of an empty slot time. Let P t be the probability 
that the channel is busy in a slot because at least one station is 
transmitting. Then, it is P t = 1 — Ilf=i(l — t^). The average 
idle slot duration can be defined as Tj = Pj ■ a = (1 — Pt) ■ cr, 



where each idle slot is assumed to have duration a. 

The average slot duration of a successful transmission. 

Consider a tagged station between the N stations in the un- 
derlined network, and let s be its index in the set {1, ... , N}. 
The probability that only the s-th tagged station is successfully 
transmitting over the channel can be defined as 
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Then, the average slot duration of a successful transmission, 
which depends on the rate-class (r) of the tagged station, can 
be evaluated as follows: 
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(14) 



+siFs + H ^y + ACK + s + DIFS 

Rc 

whereby PL is the average payload length, Hphy and Hmac 
are, respectively, the physical and MAC header sizes, r p is the 
propagation delay, DIFS is the duration of the Distributed In- 
terFrame Space, Rc is the basic data rate used for transmitting 
protocol data, and R^ is the data rate of the s-th station. 
With this setup, the average slot duration of a 

successful transmission can be evaluated as Ts = 

„(*) ( ^ u (i)\ rp(i) 



EtLiPr' Pe ) • T t 

The average collision slot duration. In a network of stations 
transmitting equal length packets with different data rates, 
the average duration Tc of a collision is largely dominated 
by the slowest transmitting stations. This phenomenon is 
called performance anomaly of 802.11b, and it has been 
firstly observed in [12]. As an example, suppose that a frame 
transmitted by a station using the rate 1 Mbps (class 1) collides 
with the packet of a station transmitting at the bit rate 11 
Mbps (class 4). Of course, both frames get lost while the 
channel appears as busy to the remaining sensing stations for 
the whole duration of the frame transmitted by the low rate 
station. Therefore, fast stations (higher classes) are penalized 
by the slow stations (low classes), causing a decrease of the 
throughput. In order to evaluate the collision probability, we 
define the class (r) collision duration as T"i r ' = H R '" Y + 
at, which takes into account the basic 
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rate Rc along with the data rate R^ of the class (r). 

For the derivations which follow, we consider a set of 
indexes which identify the stations transmitting with the r-th 
data rate: 

n(r) — {identifiers of stations belonging to rate-class (r)} 

Vr G R = {1, . . . , N R } such that J2r=i l nM l = N " I is the 
cardinality of the embraced set). With this setup, we notice 
two different collisions: 

• intra-class collisions between at least two frames belong- 
ing to the same class rate (r); 

• inter-class collisions between at least one frame of class 
(r) and at least one frame belonging to a class (j) > (r) 



As far as intra-class (r) collisions are concerned, the collision 

(r) 

probability P cl can be evaluated as follows: 
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Notice that the latter is the probability that the stations not 
belonging to the same data rate set n(r), do not transmit, times 
the probability that there are at least two stations in the same 
rate class n(r) transmitting over the channel. Notice that the 
first product within brace brackets accounts for the scenario 
in which the stations with rate in the set n(r) are silent, or 
there is only a station transmitting with rate in the set n(r). 

(r) 

As a note aside, notice that P cl = if there are no collisions 
between stations belonging to the same rate class. 

Following a similar reasoning, the inter-class (r) collision 
probability can be evaluated as: 



n d 



n 

j = r+l , 



II Ci- 
ena) 



n n a 

3 = 1 »en(j) 



which considers the scenario in which at least one station of 
class (r) and at least one station belonging to a higher rate 
class (i.e., (j) > (r)) transmit in the same slot time, while 
all the other stations belonging to lower indexed classes (i.e., 
(j) < ( r )) are silent. As a note aside, notice that pjg = if 
there are no collisions between stations belonging to different 
rate classes. 

The total class (r) collision probability is the sum of the 
previous two probabilities: 
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(17) 



while the average collision slot duration can be computed 
considering the whole set of classes r G R along with their 
collision probabilities weighted by the respective durations: 



Tc = J2 P c r) - T t ] 



(18) 



The average duration of the slot due to erroneous trans- 
missions. The average duration of the slot due to erroneous 
transmissions can be evaluated in a way similar to the one 
used for evaluating Ts and Tc'. 



N 



whereby Pj^ is defined in dT3l . and Te S> is assumed to be 

(s) 

equal to Tc since the transmitting station does not receive 
the acknowledgment before the end of the ACK timeout in 
the presence of channel errors. 
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(19) 



Average time slot duration. Given the average slot durations 
derived in the previous sections, the average duration of a slot 
time can be evaluated as follows: 



T av = T! + Tc + Ts + T E 



(20) 



TABLE I 

PHY SETUP OF THE IEEE 802. 1 IB STANDARD 



TABLE II 
Typical network parameters 



Frequency [GHz] 


"> A 
Z.4 


1 A 


A 
Z.4 


"> A 
Z.4 


Bit rate [Mbps] 


1 


2 


5.5 


1 1 


Modulation 


DBSPK 


DQPSK 


CCK 


CCK 


Chips per symbol, C s 


1 1 


1 1 


8 


8 


Bits per symbol, B s 


1 


2 


4 


8 


Channel band, B w [MHz] 


22 


22 


22 


22 


Receiver Sensitivity 


-85 


-82 


-80 


-76 


AWGN-[dBm] 











MAC header 


28 bytes 


Propag delay i~p 


— T/7s — 


PT CP Preamble 


144 bit 


PT CP Header 


48 bit 


PHY header 

1111 IJLilLIL 1 


24 bytes 


Slot timp 


20 /is 


hnQip rntf 

lul.>lL 1 dlL 


1 A/Thns 

1 ±V1.L71Jl> 


vv o 


32 


No. back-off stages, m 


5 


max 


1024 


Payload size 


1028 bytes 


SIFS 


10 /is 


ACK 


14 bytes 


DIFS 


50 fis 


ACK timeout 


364/is 


EIFS 


364 (Ms 



V. Traffic Model 

This section presents the traffic model employed in our 
setup along with the derivation of the key probabilities q 1 -^ 
and Pj S q shown in Fig. Q] The offered load related to each 
station is characterized by the parameter A^ s ^ representing the 
rate at which packets arrive at the s-th station buffer from 
the upper layers, and measured in pkt/s. The time between 
two packet arrivals is defined as interarrival time, and its 
mean value is evaluated as Jjy. One of the most commonly 
used traffic models assumes that the packet arrival process 
follows a Poisson distribution. The resulting interarrival times 
are exponentially distributed. 

In the proposed model shown in Fig. Q] we need a prob- 
ability that indicates if there is at least one packet to 
be transmitted in the queue. Probability q^ can be well 
approximated in a situation with small buffer size [7], [8], 
[9] through the following relation: 
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where T av is the expected time per slot, useful to relate the 
state of the Markov chain with the actual time spent by a 
station in each state. Such a time has been derived in (f2Qb . 
Under the hypothesis of systems employing small queues, 
probabilities q^ and Pj S q can be approximated considering 
that the probability that at least one packet arrives in the 
queue at the end of a successful packet transmission is the 
same as having at least one packet arrival in an average time 
slot duration. As a result of this simple approximation, it is 



be evaluated as follows: 
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g ( s )(a( s ) - 1) + 1 

Though simple, this approximation proved to be quite effective 
for predicting the aggregate throughput through simulation. 

VI. Physical Layer Modelling 

In a scenario with N contending stations randomly dis- 
tributed around a common AP, throughput performance de- 
pends on the channel conditions experienced by any station 
in the network. Consider a contending station at a distance d 
from the AP. Given the one-sided noise power spectral density, 
N a = —174 dBm|T=273_R" , the received SNR can be evaluated 
as [16]: 



SNR dB = P(d) |dBn 



whereby Np is the receiver noise figure (lOdB), while 
P(d) |dBm, the power received at a distance d, corresponds 
to 

P(d)\ dBm — Ptx IdBm — PL dB (24) 

Based on FCC regulations, in the 2.4GHz ISM band the trans- 
mitted power P tx |dBm amount to 20 dBm, or, equivalently, 
100 mW, while PLdB is the so-called path-loss [16]: 

PL dB = PL | dB + 10 • n p log 10 ' 



whereby PL a \ dB — — lOlogj 



G t G r \ 



The path-loss ex- 



= 10 \ (Air) 2 d'^ 

ponent, n p , depends on the specific propagation environment, 
and it ranges from 2 (free space propagation) to 3.5-4 for non- 
line-of-sight propagation, or multi-path fast fading conditions, 
in indoor environments [16]. The SNR per transmitted bit, 7, 
is defined as: 



7 | dB = SNR d 



(25) 



whereby C s stands for chips per symbol, while B s is the 
number of bits per transmitted symbol. Both C s and B s 
are summarized in Table U BER performance of the various 
transmitting modes of IEEE802.11b are shown in d26l i for 
Rayleigh fading conditions [15], [17]: 



DBPSK 
DQPSK 

CCK-5.5/11 Mbps 



2(1+7) 
1 - 
2 (°-l) 



7 



-1 (-1) 

1 



(26) 



N n - B„ 



(23) 



whereby a — 4 for 5.5 Mbps, and 8 for 1 1 Mbps, and Cf 1 = 

VII. Simulation Results and Model Validation 

This section focuses on simulation results for validating the 
theoretical models and derivations presented in the previous 
sections. We have developed a C++ simulator modelling the 
DCF protocol details in 802.11b for a specific number of 
independent transmitting stations. The simulator considers an 
Infrastructure BSS (Basic Service Set) with an AP and a 
certain number of fixed stations which communicates only 
with the AP. For the sake of simplicity, inside each station 
there are only three fundamental working levels: traffic model 
generator, MAC and PHY layers. Traffic is generated fol- 
lowing the exponential distribution for the packet interarrival 
times. Moreover, the MAC layer is managed by a state 



machine which follows the main directives specified in the 
standard [1], namely waiting times (DIFS, SIFS, EIFS), post- 
backoff, backoff, basic and RTS/CTS access mode. The typical 
MAC layer parameters for IEEE802.11b noted in Table ED [1] 
have been used for performance validation. 

For conciseness, in this paper we present a set of results 
related to following scenarios A number of 9 contending 
stations are randomly placed along a circle of radius R, while 
the AP is placed at the center of the area. Upon employing 
Equ.s (f23l>-(f25b with n p = 4 (typical of heavy faded Rayleigh 
channel conditions), we have chosen a distance R = 20 m 
in such a way that the SNR per transmitted bit is above the 
minimum sensitivity, specified in TableEl relative to the bit rate 
11 Mbps. Such stations are in saturated conditions and have 
PAR A = 8 kpkt/s. The payload size, assumed to be common 
to all the transmitting stations, is equal to 1028 bytes. In this 
scenario, another station, in the following identified as the slow 
station, is placed at 4 different distances from the AP in such a 
way that transmission occurs with the four bit rates envisaged 
within the 802.11b protocol. 

The theoretical aggregate throughput in this scenario is 
depicted in Fig. as a function of the PAR of the slow 
station. Curves in both subplots have been parameterized with 
respect to the bit rate of the slow station. Simulated points 
are noted with cross-points over the respective theoretical 
curves. The upper curves refer to ideal channel conditions, 
i.e., PER = 0, while the lower subplot represents a scenario 
in which the packets transmitted by all the stations are affected 
by a PER = 8 • 10~ 2 , which is the worst-case situation related 
to the minimum sensitivity [1]. Some considerations are in 
order. Both subplots show that the aggregate throughput is 
significantly lower than 1 1 Mbps even though all the stations 
transmit at the highest bit rate (continuous curve). Moreover, 
such a throughput reduces as a far as the PAR X s i ow increases 
reaching saturation values strongly influenced by the rate of 
the slowest station. A comparative analysis of the set of curves 
depicted in both subplots reveal the throughput reduction due 
to the presence of channel induced errors. 

VIII. Conclusions 

In this paper, we have presented a multi-dimensional Marko- 
vian state transition model characterizing the DCF behavior at 
the MAC layer of the IEEE802.il series of standards by ac- 
counting for channel induced errors and multirate transmission 
typical of fading environments, under both non-saturated and 
saturated traffic conditions. The modelling allows taking into 
consideration the impact of channel contention in throughput 
analysis which is often not considered or it is considered in a 
static mode by using a mean contention period. 

Theoretical derivations were supported by simulations. 
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